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1. INTRODUCTION 

In a global context of greenhouse gas emissions reduction, research activities are focused on the 
development of solutions, aiming at transportation-related pollution reduction. In the electric vehicles (EV) 
and mobility sector, PEMFCs are one of the promising technological advancements, for the transportation 
industry future. This technology ensures significant performances related to the EV range, compared to 
conventional architectures, powered by a hybrid energy storage system (HESS) based on SCs and batteries. 
PEMFC is an electrochemical device, used to convert dihydrogen and dioxygen to electrical power, water 
and heat. Due to its dynamic limitations, PEMFCs are activated only in the steady state. Therefore, SCs 
and/or batteries should be integrated to improve the dynamic performances and extend the PEMFC lifespan. 
This integration generates very high loads in the secondary sources during the transient regime, characterized 
by high varying energy demand. In order to ensure an optimal and more efficient contribution of fuel cell 
(FC) during transient regimes, the analysis and control of its dynamic response represents a key issue. Thus, 
several works have focused on the dynamic modeling of different energy sources, particularly FCs. These 
models take into consideration the phenomena that occur in transient regimes, such as: double layer 
capacitance, geometric capacitance and energy losses related to the system [1]-[3]. 

In this context, a simplified mathematical model of the proton exchange membrane fuel cell 
(PEMFC) has been developed, based on the variation of pressure, temperature and flow rate of hydrogen and 
air, such as the MATLAB/Simulink blocks models [4], [5]. Thus, a previous study was proposed by using the 
simplified mathematical model [6], [7], these works introduce other dynamic mechanisms, such as double 
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layer and thermal response phenomena. On the other hand, a model based on a simple approach, 
implementing a second order fuel cell transfer function is also developed in [8]. Moreover, a various studies 
have been conducted, aiming at the development of a neural network-based model, in order to predict the 
evolution of PEMFC performance [9]-[12]. 

Another dynamics models have been proposed, for example the works conducted on dynamic 
models of proton exchange membrane fuel cell (PEMFC), using electrical circuits [13]-[15]. Similarly, a 
research team [16], [17], have been particularly interested in the FC's energetic modeling, using the bond 
graph (BG) approach. The BG’s model has been evaluated using modelica software platform. On the other 
hand, a model of a PEMFC has been developed using the energetic macroscopic representation (EMR) 
approach. This work aimed at local control design, based on the EMR inversion-based rules, in order to 
control the FC’s electrical parameters [18]-[20]. 

This paper aims at the design of the PEMFC dynamic model, taking into consideration the chemical 
delay, introduced by propagation of fuel and oxidant. This work allows the construction of the FC local 
control, which integrates a chemical delay estimator. The evaluation of the developed model was carried out 
by comparing the developed model and the experimental results, obtained by BAHIA experimental platform. 
The rest of the paper is structured as follows: The second section will focus on the design of the PEMFC 
dynamic model using EMR approach. The third section is dedicated to the implementation of the FC local 
control. The last section presents the experimental results of the developed model. 


2. DYNAMIC MODEL AND LOCAL CONTROL 
2.1. PEMFC modeling 

First, a dynamic model has been developed, taking into consideration several phenomena. It’s designed 
using the EMR approach, based on causal graph principles. The developed PEMFC model is based on the 
works results of [19], [21], to which a chemical delay model has been developed and integrated, in order to 
improve the FC model behavior. This phenomenon is introduced by the propagation of fuel and oxidant [22]. 
The delay that appears during the transient regime, generates a voltage V, described in (1). Therefore, the 
implementation of the voltage V, in the Nernst equation, as shown in Figure 1, allows to improve the FC model. 
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“te.st1 

With: 

— V,: Voltage of chemical delay (V); 
— A,: Constant factor (Q) ; 

— te: Global flow delay (te = 80s). 
In the other hand, the collected equations, shown in Figure 1, describe the phenomena by domain, allowing the 
identification of the PEMFC subsystems and their interactions. Thus, the EMR model can be deduced. 
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Figure 1. Studied domains and multi-physic equations of PEMFC system 
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The EMR approach allows to establish the description and organization of multi-physical subsystems, 

based on the energy conservation principle, the different EMR blocks are collected and explained in Figure 2. 

Thus, the collected mathematical models illustrated in Figure 1, are reformulated respecting the causality rules. 

Therefore, the PEMFC EMR model illustrated in Figure 3, implements three types of blocks: 

— Source blocks: include the hydrogen source tank, ambient air, and oxygen. The cooling system and the load 
electrical sub-system; 

— Transformation blocks: in this case, the included transformation blocks describe the oxygen and hydrogen 
line losses; 

— Coupling blocks: allow multi/mono-physical coupling, i.e., fluidic, physical-chemical, thermal and electrical 
domain; 

— Accumulation block: includes the double layer capacitance, the geometric capacitance caused by the 
accumulation of the fluidic lines, and the chemical delay, induced by the propagation of oxygen and 
hydrogen, as shown in Figure 3. 
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Figure 2. Different blocs of EMR 
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Figure 3. EMR model of PEMFC 


2.2. PEMEC local control 

The objective of the local control is to ensure a system response able to track the reference, 
represented by input parameters, such as the input gas flows Qo2 ref, Qu2_ ref and the system temperature Tyc_ref 
to regulate the output voltage Vi‘ rer. Thus, the control layer is established using system's EMR model 
inversion, which allows to deduce the tuning and control chains shown in Figure 4. The adopted approach 
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leads to the organization of subsystems and the construction of the practical control structure (PCS). On the 
other hand, the research results obtained by [21], [23], dedicated on FC's modeling and control, are used in 
this study, in order to establish the local control proposed to the new FC's model. Therefore, the 
implementation of FC’s practical local control, should take into consideration the definition of the required 
estimators and reference values. 


Vecref > Va ref > En ref > AE ref —— Pc 02 ref ——— P sc 02 ref ——— 02 — Ico2 ref ——— 902 ref 


Figure 4. PEMFC control chain 


The reference parameters are defined from the equations presented in Figure 1 and illustrated in 
Table 1. The estimated parameters are calculated from the sensors as shown in Table 2. The closed-loop 
control uses a PI controller Cp; shown in (2), designed to satisfy the setpoint input represented by the oxygen 
flow demand to reach the required load voltage. In the proposed structure, the delay modeled in the model 
design section, is taken into account and estimated by the load current sensor. The Figure 5 represents the 
designed PCS based on the FC’s EMR. 
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Figure 5. PCS of the PEMFC 
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4. MATERIAL AND METHODS 

In order to evaluate the developed model, the experimental tests have been conducted on BAHIA 
bench platform as is shown in Figure 6(a). This platform is equipped with 24 PEM cells with a maximum 
power of 1.2 kW. The output voltage ranging is 13 to 27 V and nominal current is 65 A, operating in 
hydrogen/air mode [24]. The research works have been achieved on recent experimentation tests. The 
experimental platform includes three circuits: 

— The first circuit is used to assure the hydrogen supply, stored in a pressure tank; 

— The second circuit is dedicated to the oxygen supply, which is extracted from air by a controlled fan. The 
air is pumped through a membrane humidifier before being injected into FC block; 

— The third circuit is the cooling circuit, it uses water to decrease the heat of exothermic reactions, and uses 
a controller to keep the circuit regulated at a reference value of 75 °C. 

The BAHIA bench is controlled through a human machine interface (HMI) shown in Figure 6(b), 
allowing the user to introduce a power or current load profile, to set the oxygen stichometry coefficient 
ranging from 1.5 to 2.5 and the operating temperature. The measurements of the system's sensors are 
recorded in real time by an acquisition module (gas flows and pressures, temperature, voltage, current and 
power). The collected measurements are then recorded in a file in order to be processed and compared with 
the developed model response. 
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(b) 
Figure 6. Bahia bench (a) platform PEMFC and (b) human machine interface 


The validation of the proposed model should be performed using the load required power profile, in 
order to evaluate its performances. The used power supply is designed for an EV, it's calculated from urban 
driving cycle (ECE-15) and mechanical parameters of the EV in (3) [25] shown in Figure 7. The EV's 
characteristics are listed in Table 3. The power profile is scaled to 1 kW in order to be adjusted to the bench's 
maximum power supply. Therefore, the results measured on the BAHIA bench, and collected from the data 
acquisition module, are compared with the developed model results. The schematic diagram of the 
experiments is presented in Figure 8. 


M.dV 1 
Pm = V(M.g.sina+C,.M.g. cosa + —— + 5 p- S. Cy v? (3) 
Table 3. The different variables in (3) a Seay, ~<a d 
Variable Designation Value _ ES er -M 
Vv Vehicle speed (m/s) ECE-15 cycle M. g. sina ~ ’ 
a Road Grade (rad) Orad = 
Cr Coefficient of rolling 0.10 Da a 
resistance 
Cx Aerodynamic coefficient 0.30 a 
g Gravitational constant (m/s?) 9.81 m/s2 
p Air density (Kg/m3) 1.255 Kg/m3 i ; ; 7 
M Vehicle mass (Kg) 100 Kg Figure 7. Schematic of resistance forces acting on 
S Front surface m? 2.5 m2 moving vehicle 
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Figure 8. Power profile used in experimentations 


5. RESULTS AND DISCUSSIONS 

In this section, the dynamic behavior of the developed model is evaluated, comparing its results with 
those of the BAHIA bench experimental platform. The simulation model was built under 
MATLAB/SIMULINK environment, using the EMR library, developed at the University of Lille, France. 
The new PEMFC model and the local control are presented in Figure 9. In addition, the implemented load, is 
a controlled current source block, based on a programmed behavior similar to one used for the 
experimentation. The model parameter values used in this simulation are listed in Table 4. 


Figure 9. PEMFC simulation model with local control under MATLAB/SIMULINK 


The model responses are generated for a stoichiometry value of Cse 02=2 and temperature 
Ts-=75 °C. The curves shown in Figure 10, represent the responses of the developed model, for values 
defined by an empirical method of the constant factor Xe, ranging from 0 to 0.0060. Thus, according to the 
results obtained, Ae=0.004 is the value that allows to have a response approximating the response obtained 
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from the experimental bench BAHIA. In addition, according to recent research works, aiming the same 
experiments objectives, the value determined by Bahia bench, is of the same scale as the Avista Labs SR- 
12 500 -W with Ae=0.00333 Q [22] and PEMFC PC3F40 cells with Ae=0.00694 Q [26]. The constant factor 
value of Ae=0.004 Q is adopted in the developed model. 


Table 4. The different variables of the model 


Symbol Variables Value Symbol Variables Value 
F Faraday constant 96485 C/mol Rn Ohmic resistance 0.9 10° Q 
R Perfect gas constant 8.31446 J/mol/K Rt Charge Transfer Resistance 0.4 10° Q 
Po Atmospheric pressure 101325 Pa Rdey Inlet ohmic resistance for the anode 1.3 10° N.s/m? 
To Standard temperature 298.15 K Rdsy Outlet ohmic resistance for the anode 1.5 108 N.s/m> 
Ne Number of electrons exchanged 2 Cdu Fluidic capacity of the anode 4.56 10°? mN 
Ne Number of stacks 24 Rdeo Inlet ohmic resistance for the cathode 0.4 10° N.s/m? 
Sa Active Surface Area 100 cm? Rdso Outlet ohmic resistance for the cathode 1.5 10° N.s/m? 
Th Internal current O1A Cdo Fluidic capacity of the cathode 5.52 10°!° m5/N 
Ip Exchange current 4.78 10° A qu Hydrogen inlet flow 2 107 m° /s 
I Limiting current 150 A qoz Air inlet flow 5 107 m° /s 
lic Operating current 80 A he Constant factor in calculating Ed 0.0040 Q 
A Charge transfer coefficient 0.6 Te Constant time delay 80s 
900 T : 
paa [-—— Power set-point 
i [e-Test on “BAHIA” bench | 


|-æ-Simutation model A, = 0.0000 (Initial model) 
- œ -Simulation model A, = 0.0020 | | 
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Figure 10. The model responses with constant factor values Ae 


Figure 11 represents the comparison of the steady state characteristics of the proposed model, and 
the response of BAHIA bench. These characteristics were calculated using a current reference ramp function 
from 0 to 60 A. The polarization curves comparison, caused by double-layer capacitance effect and the 
undershoot phenomenon, demonstrates that the developed model agrees well with the experimental results 
during stationary regimes. 

The results shown in Figure 12, represent the PEMFC output voltage of the model and the 
experimental data, collected from Bahia bench. According to the profiles analysis, the resulting mean square 
error (MSE) of the model PEMFC output voltage, and calculated power profile, of the ECE 15 cycle is 
0.274.10°, the error found is smaller than that found in [20], which is 0.1463 shown in Figure 13. The 
comparison between the experimental data and the model developed, show that the results agree well during 
the transient and stationary regimes. 

The validation of the proposed model was also carried out by comparing the power curves of the 
proposed model and the BAHIA bench. According to Figure 14, both curves are in good agreement and have 
the same behaviour. The developed model allowed to concretize the effectiveness of the EMR approach, 
allowing to organize the mathematical models of the studied system, and the determination of the control 
parameters for Hz and O2 sources. This technique combined with the described methodology in this work, 
gave rise to a simulation platform that can be parameterized under MATLAB/SIMULINK. This simulator is 
characterized by a dynamic’s performances, which provide the possibility of its use in several domains, and 
improvement paths that can be adopted. In this case, the coupling of the PEMFC model with other sources 
such as the Battery and/or SC, constituting the (HESS). This solution allows to improve the performances of 
the main source represented by PEMFC. 
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Figure 13. Quadratic error of response to the ECE-15 cycle 
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Figure14. PEMFC power response to the ECE-15 cycle 


6. CONCLUSION 

The dynamic model proposed in this work, is in perfect agreement with the dynamic approaches of 
the FC developed using EMR. The developed model takes into consideration dynamic phenomena which 
affect its behaviour, such as fluidic interactions related to fuel flow in hydraulic circuits, the double layer and 
the charge transfer resulting from the reaction, and the voltage losses caused by electrochemical reactions in 
a cell. Furthermore, this study also integrates the local control of FC system, by using the practical inverse 
control. The developed model has been implemented and simulated in MATLAB/SIMULINK environment. 
Besides, various simulations carried out on the upgraded model, in order to validate the parameterizable 
simulation model on the BAHIA training bench, which allows to identify and characterize the control 
parameters of the FC system. The results show that the new developed model is more accurate and provides 
an excellent performance for EV design. 

Finally, the objective of this work has been focused on the modeling, simulation and validation of 
the PEMFC, which represents the main source of a HESS. Future work will focus on the design, validation 
and construction of SC and PEMFC-based HESS designed to power an EV, and the design of an energy 
management strategy ensuring the efficient and optimal supply management provided by the HESS dual- 
source. 
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